Weak topological insulator with protected gapless helical states 
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A workable model for describing dislocation lines introduced into a three-dimensional topological 
insulator is proposed. We show how fragile surface Dirac cones of a weak topological insulator evolve 
into protected gapless helical modes confined to the vicinity of dislocation line. It is demonstrated 
that surface Dirac cones of a topological insulator (either strong or weak) acquire a finite-size energy 
gap, when the surface is deformed into a cylinder penetrating the otherwise surface-less system. We 
show that when a dislocation with a non-trivial Burgers vector is introduced, the finite-size energy 
gap play the role of stabilizing the one-dimensional gapless states. 



I. INTRODUCTION 

The topological insulator has become one of the 
cutting-edge paradigms of the condensed-matter commu- 
nity since the last couple of years^— Especially high- 
lighted is the Z 2 topological insulator^ which has a band 
gap generated by spin-orbit coupling, and preserves time- 
reversal symmetry. Though the idea of Z 2 topological 
insulator stems from the two-dimensional (2D) quantum 
spin Hall effect^ its three-dimensional (3D) counterpart 
has given a stronger impact on material science, leading, 
in particular, to the reclassifying of thermo-conducting 
layered crystals such as Bi2Se3 and Bi2Te3 as "strong" 
topological insulators^ In contrast to its 2D analogue, 
the 3D Z 2 topological insulator has both weak and strong 
phases^ A strong (weak) topological insulator bears an 
odd (even) number of surface Dirac cones when it is in 
contact with the vacuum, and is characterized by a Z 2 - 
invariant Vq = 1 (vq = 0). Full characterization of a 
3D Z 2 topological insulator requires, however, a set of in 
total four Z 2 numbers: (vq;v\V2V%). 

In contrast to the topological number vq that charac- 
terizes a strong topological insulator (STI) and is asso- 
ciated with a protected surface single Dirac cone, other 
"weak indices" are generally believed to be nonrobust 
quantities. On a perfect lattice, this assertion is in- 
deed justified. A recent study, however, on the re- 
sponse of a topological insulator to the introduction of 
lattice dislocations^ e.g., screw and edge dislocations, 
suggests that such dislocation lines play the role of a 
probe for characterizing WTI, in which both strong (vo) 
and weak (V1U2 ^3) indices come into play. The authors of 
Ref4 have shown that both WTI and STI, when twisted 
by dislocations, accommodate a pair of protected one- 
dimensional (ID) helical modes. This seems to contra- 
dict the common belief that a WTI is not topologically 
robust. It is also counterintuitive that there always ap- 
pear an even number, say, two pairs of Dirac cones on 
the 2D surface of a WTI, whereas along a dislocation 
the number of protected ID Dirac cones is at most one. 
The former is susceptible to disorder, especially to inter- 
valley scattering by short-range impurities, whereas the 
latter is spin-protected from scattering by non-magnetic 




FIG. 1. A pair of screw dislocations (upper) with Burgers 
vector b — (0, 0, ±b) inserted between the two cuts (lower). 
The dislocation line is along the z-axis and parallel to the 
Burgers vector. The system is translationally invariant in the 
2-direction. 



impurities. 

The aim of this paper is to resolve the above seem- 
ingly opposing points of view on the behavior of WTI 
on a 2D surface and along a ID dislocation line. We 
propose a concrete theoretical model that is intended 
to interpolate between the two cases. To implement ei- 
ther screw or edge dislocations; see Figs. [T]and[2]), we 
first introduce two cuts extended in parallel with the z- 
axis. For analytic considerations it is more convenient 
to regard such linear cuts (of width N c ) as cylindrical 
punctures (of circumference s) penetrating the otherwise 
surfaceless system. By "cuts" we mean links on which 
electron hopping is switched off in the tight-binding de- 
scription. A pair of screw (edge) dislocations are then 
introduced around (between) these two cuts. Electrons 
in the surface states (only such electrons are relevant to 
transport characteristics) can be seen as a collection of 
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FIG. 2. Edge dislocations; a concrete implementation be- 
tween the two cuts. Here, the Burgers vector b is parallel to 
the a;-axis; b = (b, 0, 0). 

ID modes that come in pairs (Kramers' pair), moving 
up and down the punctures. These electrons also feel the 
existence of crystal dislocations. The latter plays a role 
similar to that of an (imaginary) magnetic flux piercing 
the puncture. The previously mentioned 2D and ID cases 
are naturally included within this model as the limit of, 
respectively, s — > oo and s — > 0. We follow the evo- 
lution of electronic states along such punctures with a 
non-trivial lattice distortion as s is varied. It is revealed 
that the topological stability of protected ID gapless he- 
lical modes stems from a finite-size energy gap associated 
with the spin Berry phase. The latter has been a subject 
of much theoretical attention^r— in the context of pecu- 
liar Aharonov-Bohm oscillations observed recently in a 
system of STlJ^ The protected gapless modes along dis- 
location lines have been studied also from the viewpoint 
of engineering thermo-electric materials^ 

II. MODEL 

In the bulk (outside the punctures and away from the 
dislocation) we consider a lattice version of the following 
simplified model for 3D Z 2 topological insulator 

H = Ak^ + (A - Sfc^ho (1) 

where repeated indices should be summed over /j = 
1,2,3. 7-matrices are chosen, e.g., as, 

7m = T z^n, 7o = t x (2) 

for /j = 1,2,3. Then, following the same type of pro- 
cedure as described in Refs . 14 ' 15 , we place the model on 
a 3D square lattice of size N x x N y x N z , and impose, 
unless stated otherwise a periodic boundary condition in 
each direction. 



Away from the two cuts and dislocations, our tight- 
binding Hamiltonian reads, 

H = ^ {m(k z )\x,y,k z )(x,y,k z \ 

x,y,k z 

+ (t x \x + l,y,k z )(x,y,k z \ 

+ t y \x,y+ l,k z )(x,y, k z \ +h.c.)} (3) 

where 

m(k z ) = Asmk z r z a z + (A — 6B + 2B cos k z )r x , 
A 

tx = l-^T z (T x + Bt x , 

l y = i -2 Tz(T y + BTx ' ( 4 ) 



A. Cuts 

In order to implement a punctured geometry and in- 
troduce dislocations on the square lattice, we first de- 
form the punctures into the form of a "cut" (see the 
lower panel of FIG. [T) of length N c (its circumference 
is s — 2N C ). We introduce two cuts, then a pair of screw 
(Fig. [T|) or edge (Fig. [5]) dislocations between them. As 
shown in these figures, here the two cuts are placed along 
the z-axis, and between the two crystal layers: y = and 
y = 1 (as well as between y = N y /2 — 1 and y = N y /2). 
Between these crystal layers hopping is turned off for 
x = 1, ■ ■ ■ , N c . Introduction of these two cuts breaks the 
discrete translational invariance (crystal periodicity) in 
the (x, y)-plane, whereas it preserves the translational in- 
variance in the z-direction, i.e., k z is still a good quantum 
number. In the following, we will extensively investigate 
energy spectra: E = E(k z ) of the system in the presence 
of screw or edge dislocations. 



B. Screw vs. edge dislocations 

1. Case of screw dislocations 

A pair of screw dislocations can be introduced between 
the two cuts by dislocating the hopping matrix elements 
in the region between the two cuts (Figs. [TJ, i.e., for y = 
1, • • • ,N y /2, say, between the two crystal layers x = 
and x = 1 as 

t x \x + l,y,z)(x,y,z\ t x \x + l,y,z - b){x,y,z\. (5) 

b measures the strength of the dislocation, i.e., the mag- 
nitude of the Burgers vector. This is equivalent to 
" twisting" the same hopping matrix elements by a factor 
e ik z b m £ ne £ 2 _diagonalized basis.— Note that the cuts 
and twist structure is translationally invariant in the z- 
dircction, and k z is still a good quantum number. 
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TABLE I. Three distinct topological phases of the square lat- 
tice Dirac mode: its low-energy effective Hamiltonian around 
the F-point is given in Eqs. (HJ , (J2j) . 



So Si S2 83 (vq- V\ViVz) phase 

A/B < + + + + (0;000) trivial 

0<A/B<4 - + + + (f;000) STf 

4 < A/B < 8 + + (0; 111) WTI 

8 < A/B < 12 + (1; 111) STI 

12 < A/B (0;000) trivial 



2. Case of edge dislocations 

In the case of (a pair of) edge dislocations with burg- 
ers vector b = (±&, 0,0), we suppress b nearest-neighbor 
hopping amplitudes in the x-direction between x = xq 
and x — xo + 1 + b and also between the two cuts 
(1 < V < Ny/2), and instead introduce a "skipping" 
process, 

t x \x + 1 + b,y,z)(x ,y,z\, (6) 
to the tight-binding Hamiltonian 



C. Strong vs. weak indices 

The 3D topological insulator model we employ has 
three distinct topological phases as shown in Table 1. 
Indices So, S\, S2, S3 in the Table are parity eigenvalues: 
±1, respectively, at T = (0,0,0) (for So), at three incquiv- 
alent but symmetric Appoints: (it, 0, 0), (0, it, 0), (0, 0, tt) 
(for Si), at three M-points: (0,7r, tt), (tt, 0,7t), (tt, tt,0) 
(for S2) and at R = (tt,tt,tt) (for S3). The above eight 
(by distinguishing incquivalent points) symmetry points 
(r, X, M and R) are also called time- reversal invariant 
momenta (TRIM) of the 3D Brillouin zone. These parity 
eigenvalues arc in our model related to the strong and 
weak indices as, 

(-I)" = S S 1 S[S'{S 2 S 2 5 2 % = S SlS 3 2 S 3 , (7) 

(-1)" 1 = W 2 % = 5i5l5 3 , (8) 

= S[S 2 S' 2 % = 3,3%, (9) 

(_!)«. = S'IS 2 5' 2 S 3 = S 1 5lS 3 . (10) 

Here, we have distinguished, for later convenience, three 
5 2 's and S 3 s at symmetric but inequivalcnt TRIM (the 
value of these <5's are identical in our model with high 
symmetry; as for definitions of these <5's, see Fig. [8]). 

In the following, we focus on the WTI phase: 4 < 
A/B < 8, and study how a protected ID helical pair 
arises from a topologically fragile surface of a WTI. 




FIG. 3. Calculated spectrum of surface Dirac cones, in the 
WTI phase (A/B = 6, A = B = 1). Periodic boundary 
condition in the ^-direction is relaxed; the system forms a 
slab, or a torus of finite thickness. The two Dirac cones are 
located at TRIM: (tt,0) and (0,vr). 



III. ENERGY SPECTRUM OF WTI IN THE 
PRESENCE OF PUNCTURES AND 
DISLOCATION LINES 

A WTI has an even number of Dirac cones on its sur- 
face as depicted in Fig. [3l Here, the surface is chosen 
normal to the j/-axis, i.e., a WTI occupying the half space 
y < is in contact with the vacuum occupying the re- 
maining half at the y = surface. The two Dirac cones 
are located at two TRIM's: (0,7r) and (tt, 0) in the sur- 
face coordinates (k z ,k x ). 

A. Finite-size energy gap of surface Dirac cones on 
a cylindrical surface 

Imagine deforming this flat surface into a cylindrical 
tube. The tube is further deformed adiabatically into 
a cut of Fig. Q] and Fig. [5J The two Dirac cones are 
now projected onto the fc 2 -axis, as shown in the Fig. [4] 
Notice that the two projected Dirac cones at k z = 
and k z = tt have acquired a finite size gap in the up- 
per panel. Note that here the twist is not introduced yet. 
The appearance of a gap is a rather unexpected phe- 
nomenon, if one recalls that carbon nanotubes become 
either metallic or semiconducting depending on the way 
a graphene is rolled up into a tubeJ^ Here, a crucial dif- 
ference from the carbon nanotube case is that the Dirac 
cone involves a real spin and not a sub-lattice pseudo- 
spin. The procedure of rolling up a fiat surface into a 
tube introduces a 27r-rotation in spin space along a con- 
tour winding around the tube once^— The resulting —1 
factor changes the boundary condition around the tube 
from periodic to anti-periodic: 

e *( P *+fc< 0) ) s x = 1 (11) 

Here, we have decomposed the total crystal momentum of 
an electron into short- and long-wavelength components: 

k = k {0) +p. (12) 
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FIG. 4. Energy spectrum E(k) of WTI (A/B = 6) in the 
presence of screw dislocations. Here, the ID momentum k is 
chosen to be along the cuts; k — k z . The Burgers vectors 
are gives as b = (0, 0, ±6) with b = (dislocation is absent), 
b = 1, b = 2 and 6 = 3 respectively, in the top, second, third 
and bottom panels. The calculation is done for a system of 
size, N x x N y = 16 x 16, and cut width iV c = 8. The other 
parameters are set as A — B — 1. 



p = ipxiPy) refers to the long- wavelength component 
measured from the Dirac point. The short- wavelength 

component = (ki°\ky ^) is, on the other hand, 

a crystal momentum at the Dirac point, and typically 
= 0, 7T. Recall here that the circumference s = 2N C 
of the cut is, by its construction, an even integer multiple 
of the lattice constant, since the cut is made by discon- 
necting N c links of an otherwise locally perfect crystal. 
This signifies that 

e* k * )s = 1 (13) 

always holds. As a result, the anti-periodicity of the 
boundary condition (fTTj) must be taken care of solely 
by the long-wavelength part of the crystal momentum, 
and eliminates, as seen e.g. in the spectrum of Fig. [4] 
(top panel), states on the line p x = crossing the very 
bottom of a Dirac cone. Low-energy states in the same 
figure consist of p x = ±7r/s, leading to occurrence of a 
finite-size gap, 

AE = 2A-. (14) 
s 

in the spectrum. 



B. Screw dislocations 

The second panel of Fig. [4] shows, on the other hand, 
the spectrum when the system is twisted by a pair of 
screw dislocations with Burgers vector b = (0, 0, ±6) 
where 6=1. Such a lattice scale deformation modi- 
fies the periodicity of the wave function associated with 
the short-wavelength component of the crystal momen- 
tum, i.e., ki°^ = 7r in the present case. Note that the 
entire effect of a screw dislocation can be concentrated 
on hopping amplitudes across a single surface, as in Eq. 
([5]). Its influences on the electronic wave function sums 
up to a phase shift e lkzb on crossing the same surface 
(here, this is a surface inserted between the two crystal 
layers y = and y = 1). Thus, adding this phase shift 
to (fTTj) and taking Eq. (fl~3|) into account, one finds that 
the appropriate boundary condition in the presence of a 
screw dislocation reads, 

e i P *s x e ik?h x ^ = L ( 15 ) 

Note that here a small additional phase factor e lPzb which 
modifies only gapped solutions with p z ^ has been omit- 
ted for the sake of clarity. Eq. (|15p dictates that only 
the surface Dirac cone projected onto = tt is sus- 
ceptible to the change of the magnitude of the Burgers 
vector, and closes the gap (i.e. the p x = state is now 
allowed^) when b is an odd integer. 

Some examples confirming this even/odd feature are 
shown in Fig. 0J In the last two panels of Fig. [U one 
can also observe that Kramers pairs at k z = exchange 
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FIG. 5. Similar plots as the last two panels of Fig. [4] here 
in the case of edge dislocations. The Burgers vector is here 
chosen as b — (±6, 0, 0) with b = 2 and b = 3, respectively, in 
the upper and lower panels. 
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FIG. 6. Protected gapless helical modes along a pair of screw 
dislocations (b = 1). The spectrum is calculated at A/B = 6 
for two different values of the cut- width N c ; N c = (upper) 
and N c = 16 (lower). The size of the system is chosen as 
N x x N y — 32 x 8. See the text for discussions. 



their partners as k z evolves up to k z = w, in accordance 
with the twisting of boundary condition. 

The anti-periodic boundary condition the re- 

sulting finite-size gap (fT4)) , as well as the twisting of 
the boundary condition such as Eq. (fT5|) also under- 
lie the origin of the anomalous Aharonov-Bohm oscilla- 
tions observed recently in Bi2Se3 nanoribbbons.— In the 
Aharonov-Bohm geometry, the twisting of the bound- 
ary condition a la Eq. (fl5|) is caused, not by a dislo- 
cation, but instead by a magnetic flux tube penetrating 
the punctured— 



C. Edge dislocations 

The above argument needs to be modified in the case 
of an edge dislocation associated with the same dislo- 
cation line. Such defects can be introduced e.g., as in 
Fig. [2j in which dislocations terminate at a cut of fi- 
nite width similarly to the case of a screw dislocation. 
The Burgers vector in this implementation is along the 
x-axis, b = (b, 0, 0). Here, b is the number of subtracted 
(added) layers between the two cuts. Recall that for an 
edge (screw) dislocation the Burgers vector is perpen- 



dicular (parallel) to the dislocation line (parallel to the 
z-axis here) . The effect of such an edge dislocation on the 
electronic wave function can be fully taken into account 
as a change of the boundary condition (jTTJ) , i.e., by the 
replacement: s — > s + b in the same equation. This leads, 
when (|13p is taken into account, to 

e lp ' s x e iki ° )b x (-1) = 1, (16) 

i.e., a twisted boundary condition analogous to Eq. (|15[) 
but with k z ^ replaced by k x °\ Note that again there is 
a small additional phase factor e* Pxh which appears only 
when p x 7^ 0. Eq. (fTo) dictates, in contrast to Eq. (fT5]) . 
that among the two surface Dirac cones projected onto 
fcz-axis, only the one with p x — n is susceptible to the 
presence of edge dislocation, and closes its finite-size gap 
when b is an odd integer. Two panels of Fig. [5] indeed 
confirm this even/odd feature in a few non-trivial cases: 
6 = 2,3. Notice that protected gapless modes appear at 
k z = [6 = 3 (odd) case] , in contrast to the case of screw 
dislocation. This is because here the underlying surface 
Dirac cone responsible for the gap closing is located at 
(k z ,k x ) = (0,7r), projected naturally to k z = 0. 
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relevant for ^ — 2. 
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FIG. 8. Parity eigenvalues determining both the indices 
(uo; U1U2V3) and the position of surface Dirac cones on the 
projected 2D-plane; here chosen as the (k z , fc^-plane. Case 
1: screw dislocation [6 = (0, 0, 6)] — 5" 82 and S2S3 occurring 
at (k z , k x ) = (n, 0) and (fc z , k x ) — (ir, n) are relevant. Case 2: 
edge dislocation [6 — (6,0,0)] — S1S2 and 8263 occurring at 
(k z , k x ) = (0, 7r) and (fe z , k x ) = (n, n) are relevant. 



FIG. 7. Protected gapless helical modes along a pair of edge 
dislocations (6 = 1). The spectrum is calculated at A/B = 6 
for two different values of the cut- width N c ; N c = 1 (upper) 
and No = 16 (lower). N x x N y = 32 x 8. 



IV. FINITE SIZE GAP OF PROJECTED 2D 
DIRAC CONES AND THE PROTECTED ID 
GAPLESS HELICAL MODES 

We have seen that 2D surface Dirac cones attains a 
finite-size mass gap, when the surface is deformed into 
a tube of finite circumference s (cf. Eqs. (|TT|) . (fT4l) ). 
We point out that this observation is the key for under- 
standing the mechanism of how the originally fragile 2D 
surface Dirac cones of WTI acquires robustness upon the 
introduction of a dislocation, and transforms into pro- 
tected ID gapless helical modes. 

In the absence or presence of trivial (b: even) dislo- 
cations, the finite-size gap evolves continuously into the 
bulk gap as s — > 0. When b is odd, the same evolution 
gives robustness to the gapless modes. When the cir- 
cumference s of the puncture, around which the crystal 
is dislocated, is finite, the gapless modes are separated 
from the (gapped) continuum only by an energy of order 
A/s. As the size of the puncture is reduced, only the gap- 
less pairs stay intact, and its unique property that it is 
topologically protected manifests, making it distinguish- 
able from the rest of the spectrum. Projected Dirac cones 
without a pair of protected ID gapless helical modes be- 



come indistinguishable from the gapped bulk spectrum. 

Fig. |5] (Fig. [7]) depicts such an evolution in the pres- 
ence of screw (edge) dislocations. In the two figures, one 
can observe, upon reducing the size of the cuts (from 
lower to upper panels) from N c = 16 either to N c = 
(screw case) or to N c = 1 (edge case), that the gapless 
helical pair isolates. Note that in the case of edge dis- 
locations one cannot reduce the cut width smaller than 
b. Note also that in these plots separation between the 
two cuts is relatively small (N y /2 = 4) in order to take 
the width of the cut sufficiently large (jV c = 16). For 
this reason, there appears a finite interference between 
the ideally gapless counter-propagating modes, each lo- 
calized in the vicinity of two dislocation lines. Of course, 
when the size of the cut is finite, the wave function of 
the gapless mode is extended almost uniformly around 
the cut. The wave function shows a sharp peak in its 
amplitude in the vicinity of a dislocation line in the limit 
N c — > (not shown). 



V. RELATION BETWEEN WEAK INDICES 
AND PROTECTED ID HELICAL MODES 

What is the relation between the weak indices and con- 
dition for the appearance of protected ID helical modes? 
A deep connection between these two a priori unrelated 
quantities becomes manifest when expressing both the 
weak indices and the latter condition in terms of the par- 
ity eigenvalues at the eight bulk TRIM's, since our sys- 
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FIG. 9. (a) Necessary arrangement of the 2D indices intro- 
duced in Fig. [8] for the appearance of protected ID helical 
modes along a screw dislocation in the z-direction. Only rel- 
ative signs are relevant. Column (b): WTI example satisfying 
the condition in (a); 4 < A/B < 8. Column (c): STI example 
satisfying the condition in (a); 8 < A/B < 12. 



tern has inversion symmetry^ The expressions for weak 
indices in terms of So , 5\, 82 and £3 were given in Eq. ()10[) . 
In order to identify the condition for the appearance of 
protected ID helical modes in terms of So, Si, 82, 83, and 
<5's at symmetric but inequivalent points, we project the 
3D reciprocal space in which the eight 3D TRIM are de- 
fined in two steps; first onto the 2D reciprocal surface on 
which surface Dirac cones appear, then further onto ID 
/cz-axis on which protected ID helical modes appear. Fig. 
[8] shows how the eight parity eigenvalues (among which 
only four are independent) determine the weak indices, 
say, V3, upon projected onto the (k z , fc^-plane. Products 
of two indices at four 2D TRIM determine the position 
where surface Dirac cones appear^ 

Fig. [5] shows, on the other hand, that the appearance 
or disappearance of protected ID helical modes along a 
screw dislocation in the z-direction is related to a relative 
sign of indices, 8" 62 and S'283 occurring at (k z ,k x ) = 
(■7T, 0) and (k z ,k x ) = (ir,ir). When these indices have 
opposite signs, 



(S'{5 2 ) x (S' 2 S 3 ) 



-1, 



(17) 



face Dirac cone that is projected to k z — tt. This pro- 
jected Dirac cone acquires a finite-size mass gap that is 
susceptible to the change of the boundary condition (cf. 
Eq. (fT5]0 caused by the twisting associated with, e.g., 
a screw dislocation. The projected Dirac valley features 
a protected ID helical modes when b is an odd integer. 
Notice, on the other hand, that the same combination 
of parity eigenvalues as Eq. (|17[) has appeared in Eq. 
pOll (see also Fig. [8j). Thus, "1/3 = 1 and b is odd" - 
(A), is both a necessary and sufficient condition for the 
appearance of protected ID gapless helical modes. 

The situation is different for an edge dislocation, where 
the dislocation line is taken to be parallel to the z-axis 
but with a Burgers vector b = (±6,0,0). In this case, 
the appearance of protected ID modes is related to a 
relative sign of the indices 8182 and 8' 2 83 occurring at 
(k z ,k x ) = (0,7r) and (k z ,k x ) — (tt,tt); see Fig. [8j When 
these indices have opposite signs, i.e., 



(S^) x (S' 2 S 3 ) 



-1, 



(18) 



an odd number of surface Dirac cones are susceptible to 
the change of boundary condition (flU)) associated with 
the insertion or subtraction of crystal layers between the 
two cuts. The same combination of S's as Eq. psp has 
appeared, in contrast to the previous case, in Eq. (J8]) (see 
also Fig. [H]). Thus protected ID gapless modes appear 
in the present case iff "v\ = 1 and b is odd" — (B). 

The above statements (A) and (B) concerning the ap- 
pearance of protected ID gapless modes are consistent 
with the expression, 

M-b = ir mod 27r, (19) 

which has appeared in Ref4. This formula is a straight- 
forward generalization of the criteria (A) and (B) to the 
case of the absence of inversion symmetry. In Eq. CG 
the vector M is defined as, 



M = + V2G2 + ^ 3 G 3 ), 



(20) 



there appears an odd number of, i.e., a single 2D sur- 



in terms of reciprocal lattice vectors, G\, G2 and G3. 
Note that the same formula can be derived by considering 
winding properties of a Bloch electron in an extended 
parameter space incorporating the dislocation linesJ^ 

In Fig. [9](a) the lower-left index S081 is irrelevant, since 
Dirac cones projected onto k z = is insensitive to the 
change of boundary condition (cf. Eq. (fll?]0 . This means 
that the dislocation probes only weak indices. Protected 
ID gapless helical modes similarly appear both in the 
WTI and STI phases with the same weak indices (see 
Fig. [S] columns (b-c)). 

Fig. [10] shows STI example on the (dis) appearance 
of protected gapless modes, both in the screw and edge 
dislocation cases. Recall that in the WTI phase protected 
gapless modes along an edge dislocation appear at k z = 0, 
whereas here the same protected modes appear at k z = ir, 
even though the two phases are characterized by the same 
weak indices; {v ;v 1 v 2 v 3 ) = (0; 111) [WTI: 4 < A/B < 
8] and (vo;viV2Vz) = (1; 111) [STI: 8 < A/B < 12]. 
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FIG. 10. Plots similar to Fig. H and Fig. [5] here in the STI 
phase {A/B = 10). The upper (lower) panel corresponds to 
the case of screw (edge) dislocations, b is chosen as, b = 3 and 
b — 2, respectively, in the two cases. 



VI. CONCLUSIONS 



In this paper, we have addressed the question: how 
"weak" is a WTI? The existence of protected gapless he- 
lical states parasitic to a dislocation line of a WTI seems 
per se contradictory to the fragility of the even numbers 
of surface Dirac cones of a WTI. Using a simple model 
for a topological insulator implemented on a square lat- 
tice, we have systematically studied the nature of elec- 
tronic states in the presence of dislocation lines. In order 
to resolve the apparent contradiction between the stabil- 
ity of ID gapless helical modes and the nonrobustness 
of 2D surface Dirac cones, we have invented and stud- 
ied a modified variant of the defect-free model in which 
a dislocation is extended along a cylinder of finite cir- 
cumference. The unexpected stability of the ID helical 
states was identified as an interplay of the finite-size en- 
ergy gap specific to surface states of a 3D topological 
insulator and twisting of the boundary condition due to 
topologically nontrivial geometry. This scenario is closely 
related to the mechanism of recently observed anomalous 
Aharonov-Bohm oscillations in a STI of ribbon geometry. 
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